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The yield of amorphous solids like metallic glasses under external stress was discussed asserting 
that it is related to the glass transition by increasing temperature, or that it can be understood using 
statistical theories of various sorts. Here we study the approach to stress-controlled yield and argue 
that neither assertions can be supported, at least at low temperatures. The yield of amorphous 
solids at low temperatures is a highly structured phenomenon, characterized by a specific series of 
mechanical instabilities, and having no similarity at all to fluidization by increased temperature, real 
or Active. The series of instabilities followed by stress controlled yield at low but finite temperature 
protocols can be predicted by analyzing athermal quasi-static strain controlled protocols, making 
the latter highly relevant for the deep understanding of the mechanical properties of amorphous 
solids. 


The lack of crystalline order in amorphous solids makes 
the study of their mechanical properties a very challeng¬ 
ing subject of high current interest. One particularly 
intriguing characteristic is their yield-stress ay which is 
defined as the lowest external mechanical stress above 
which unbounded plastic flow can be reached M- For 
lower values of the external stress the amorphous solid 
deforms but remains able to develop an internal stress 
that balances the external stress, leaving the system in 
mechanical equilibrium. Naturally one wants to under¬ 
stand the yielding process, sometime with the hope of 
learning how to improve the material atomic composition 
or its protocol of preparation. In this Letter we address 
the atomistic theory of the approach to yield under stress 
control protocols. 

Scanning the literature on the dynamics of yielding 
of amorphous solids one finds basically two schools of 
thought. One school attempts to describe yielding using 
one or another sort of statistical theory that assumes that 
the process can be modeled as a random series of creation 
and annihilation of soft spots, or regions of plasticity that 
are sometime referred to as shear transformation zones 
g. Another school of thought asserts that mechanical 
yielding under stress can be understood as a kind of flu¬ 
idization, similar to the inverse glass transition due to 
the increase of temperature g. Both schools of thought 
may or may not invoke a notion of Active or effective tem¬ 
perature which, in order to explain the yield, is assumed 
to increase even though the thermodynamic temperature 
remains low g. We will provide strong evidence in this 
Letter that yielding is a highly structured, not random 
at all, series of instabilities that can be predicted a-priori 
from zero-temperature considerations. Moreover, we will 
show that that there is no increase in Active temperature 
near yield, and there is no necessity to invoke such an 
increase. The phenomenon of yield is purely mechani¬ 
cal and is easy to understand once the right theoretical 
framework is provided. 

To understand the phenomenon of stress controlled 
yield one requires a good method to observe the phe¬ 
nomenon with a high degree of accuracy. In recent work 


(where full details can be found) such a simulational 
method was described using a “variable shape Monte- 
Carlo” technique g. This method introduces strain into 
the simulation box by Arst deAning a square box of unit 
area where the particles are at positions Next one 
deAnes a linear transformation h , taking the particles to 
positions r; via r j = h • S;. In order to prevent rotations 
of the simulation box, the matrix h should be symmetric. 
The current area of a system becomes the determinant 
V =| h |. To proceed, two kinds of trial moves are con¬ 
sidered: in the Arst we perform n standard Monte Carlo 
moves (displacement of the particle positions given by s)) 

s new = s old + Ss , l <i < N. (1) 

In this equation the a component of the displacement 
vector of a particle is given by 

Ss a = A W (2£“ - 1), (2) 

where A s max is the maximum displacement and £“ is an 
independent random number uniformly distributed be¬ 
tween 0 and 1. After n sweeps deAned by Eq. m the 
transformation h changes according to 

h new = h old + Sh , (3) 

where elements of the random symmetric matrix Sh are 
deAned by 


Shij — Ah max (2^ij 1), i ^ j. (4) 


Here A h max is the maximum allowed change of a ma¬ 
trix element and is an independent random number 
uniformly distributed between 0 and 1. The value of 
A h ma x and the maximum displacement of particle posi¬ 
tions A s max are selected to obtain a desired acceptance 
rate (usually between 30% and 50%). For each kind of 
move the trial conAguration is accepted with probability 


Ptr = min 


1, exp 



( 5 ) 
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arrests it again reaches detailed balance. (Sometime dur¬ 
ing the equilibration process the system resides in config¬ 
urations that are part of a metastable state long enough 
to reach “almost-detailed-balance” also there, before dis¬ 
covering the escape channel). Panel (d) shows that for 
a xy * =0.5 the system collapses mechanically and distorts 
indefinitely with 7 increasing without any obvious limit. 
Thus we can conclude that for the present realization at 
a temperature of T = 0.005 the yield stress lies between 
0.46 < cry < 0.5. 

At this point we can ask whether approaching the yield 
point from below is similar in any way to approaching the 
glass transition point by increasing the temperature. To 
this aim we study the Hessian matrix which is defined as 


Ra ,p = d 2 U(r 1 ,r 2 ,---r N ) 
drfdr f 

% J 


(9) 


FIG. 1: Internal stress a xy vs. strain 7 in one realization 
for four different values of external stress at a tem¬ 

perature T = 0.005 in Lennard Jones units. Panel a-d: 
er^y* = 0.44, 0.45, 0.46, 0.5 respectively. 


Here AG is the enthalpy change due to the move, 

AG' = U (7 + <57, rT w ) - U( 7 , r° ld ) - Va e x f8 1 . ( 6 ) 

where 7 is the system’s strain. Temperature is mea¬ 
sured here in Lennard-Jones units with the Boltzmann 
constant taken as unity. To specialize the technique to 
stress-controlled simple shear simulations at any desired 
temperature one chooses the following h matrix 

h =r(;i). u) 


where L is the length of the square simulation box and 
7 is the simple shear strain, with the volume of the sys¬ 
tem V=L 2 being conserved. The external stress in this 
protocol is given by 


_ ext 

(T 


0 ^ x f\ 

< X y 0 ) 


( 8 ) 


The nature of the results of the Monte Carlo protocol 
is presented in Fig [1] In panel (a) the external stress is 
a xy t = 9-44 and the system equilibrates with stress fluc¬ 
tuations around the value of the external stress. The red 
cloud of points represents the values of 7 and a xy of the 
configurations that are visited during the Monte Carlo 
process. The system reaches detailed balance under the 
action of the enthalpy [ Eq. ©]■ Panels (b) and (c) show 
that for slightly higher external stress = 0.45 and 
a xy = 9-46 the system visits several metastable states 
before arresting again at values of the internal stress that 
fluctuate around the external stress. When the system 


At mechanical equilibrium and T = 0, The Hessian ma¬ 
trix has 2 N eigenvalues that are all positive. At finite 
temperature one needs to arrest a given configuration 
and compute the Hessian H for the arrested positions. 
Then the spectrum contains negative eigenvalues whose 
number n neg increases with the temperature 00 . If 
the approach to yield were analogous to the increase in 
temperature towards the glass transition, one would ex¬ 
pect n neg to increase as we approach a yield point. In 
our simulations we can arrest the Monte Carlo protocol 
at any value of < 7 ®“* and any value of T for which the 
system reaches detailed balance with (<r xy ) = cr®^ 4 . We 
then compute the Hessian matrix for these arrested con¬ 
figurations with the results shown in Fig. [2] We find 
that the ratio of the number of negative eigenvalues to 
2TV is fairly independent of N and we thus report the 
number r neg = n neg /2N. In the upper panel we present 
the average over 500 realizations of the ratio of negative 
eigenvalues r neg for zero external stress (black line) and 
for various external stresses up to the yield stress as a 
function of temperature. One sees that r neg does not de¬ 
pend at all on the external stress. The system does not 
“prepare itself’ for failure by increasing a Active temper¬ 
ature. In the lower panel the same conclusion is drawn 
from the independence of r neg on the external stress, here 
shown for different temperatures. 

It remains therefore to understand what is then the 
phenomenon of yield. It is very helpful to consider, along¬ 
side the stress-controlled Monte Carlo simulations, also 
an athermal quasistatic strain-control protocol in which 
the strain is increased in steps, forcing the system to re¬ 
turn all the time to mechanical equilibrium at T = 0. An 
example of the stress-strain dependence in such a “AQS” 
protocol is shown in panel (a) of Fig. [3] The elastic in¬ 
creases in stress are punctuated by plastic instabilities 
with sharp decreases in stress (and energy). After ev¬ 
ery plastic instability the system is forced to return to 
mechanical equilibrium using energy minimization. Thus 
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FIG. 2: Average over 500 realizations of the dependence of 
the fraction of negative eigenvalues of Hessian matrix r ne g as 
a function of stress and temperature. Upper panel: different 
external stresses as a function of temperature. Lower panel: 
different temperatures as a function of external stress. 




FIG. 3: Panel (a) : stress vs. strain in an AQS protocol that 
forces the system to follow the lowest energy path. Panel (b): 
the Monte Carlo protocol with = 0.45 and T = 0.005. 
Panel (c): the Monte Carlo protocol with cr™ 4 = 0.46 and 
T = 0.005. 


the AQS protocol describes a deformation along the min¬ 
imal energy path of a very complex energy landscape that 
is being distorted by the external strain. The system fol¬ 
lows the series of saddle node bifurcations that are caused 
by this distortion keeping itself always at the lowest pos¬ 
sible energy. The saddle bifurcations result in shifting 
the system from one local minimum of the landscape to 
another. Panels (b) and (c) show that the Monte Carlo 
stress-controlled protocol at finite temperature follows es¬ 
sentially the same path. Imposing an external stress the 
system attempts to balance it by increasing its internal 
stress, making stops in the metastable stress-strain val¬ 
ues that correspond to the deepest basins of the energy 
landscape which are exposed by the AQS plastic drops. 
In other words, the process of “yield” is neither statisti¬ 
cal nor random, it follows carefully the easiest path on 
the energy landscape. 

This picture indicates that in fact there is no single 
“yield-stress” but rather a series of such, and we can or¬ 
der them as <J^\ etc. in order of increasing heights 
of the stress peaks in the first panel of Fig[3] Of course, 
temperature fluctuations can cause “glitches” in which 
a the system may choose a sub-leading local minimum 
from which the path may deviate from the AQS analog. 
Nevertheless the finite temperature path is crucially sen¬ 
sitive to the topography of the energy landscape greater 
basins. 

As additional demonstration that the stress-controlled 
yield is dominated by the underlying topography we turn 
now to the temperature dependence of the observed yield 
stress 0 -O) (J 1 ). We expect that if the yield follows the 
same saddle-node instabilities that are so well known 
from AQS protocols, this must be reflected in the scaling 
law of A a[ n) (T) = ct£*>(T = 0) - The scal¬ 

ing law is easily predicted based on our knowledge that 
during a saddle node bifurcation as 7 —> yp the criti¬ 
cal eigenvalue Ap of the Hessian matrix goes to zero like 
V'yp ~ 7 HH3- To find the height of the energy barrier 
we can expand the energy in the direction of the softening 
mode to obtain 

U(S) = U miD + ^Ap( 7 )S 2 + \l<s 3 + .... (10) 

The reducing maximum is located at value of S for which 
dU / dS = 0 or S ~ y/Xp. From here we find immediately 
that the energy barrier reduces like A U ~ (y/X p ) 3 ~ 
[ 7 p — y] 3 / 2 . At T = 0 the barrier cannot be overcome 
until A U = 0. But for finite T we expect that yield can 
occur when T ~ A U. We thus conclude that we expect 
a scaling law (taking Act ~ Ay): 

A a ( f\T) = CT 2/z . (11) 

This law is expected to hold every time that the exter¬ 
nal stress in the Monte Carlo stress-control protocol ap¬ 
proaches a value of ct^, which can be associated with 
some higher peak in the stress vs. strain graph of the 
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FIG. 4: Temperature dependence of A a^(T) for the first 
yield (corresponding to the first peak in the AQS strain vs. 
stress curve), for two different realizations (squares/triangles). 
The slope in this log-log plot agrees well with the scaling 
exponents 2/3 of Eq. (Hill . 

AQS strain control protocol. Fig. [3] panel (a). Indeed, 
in Fig. [4] we compare the scaling law (fill) to the mea¬ 
sured yield stress for different temperatures associated 
with overcoming the first peak, i.e. for A a^(T) . Of 
course, the same scaling law is expected to occur close to 
yield for all the critical values a ^, including the highest 
which is the classically considered “yield stress”. 

In general, the number and structure of those stress 
drops and energy landscape basins corresponding to 
cr^.-.a^, depends on the glass preparation protocol, 
e.g. quench rate, and on system size. For slow quench 
rates and for bigger system one can find that the first 
stress-drop also corresponds to the highest peak, and 
to the final collapse of the sample, and thus uniting the 
two definitions for yielding. 


This 2/3 scaling law is also similar to the earlier re¬ 
ported strain-controlled El and experimental results 
|12j |. supporting the mechanical nature of this yielding 
phenomenon and its relation to the saddle-node bifurca¬ 
tions of the energy landscape. 

In summary, we have shown here that stress-controlled 
yield in amorphous solids is akin to strain-controlled AQS 
protocols in the sense that the yielding process follows 
(at low enough temperatures) the minimal path to in¬ 
creased strain, along the “bottom” of the energy land¬ 
scape seeking the lowest saddle nodes that open up for 
increased strain. The AQS analog has no “yield” since 
one forces the system to return to mechanical equilib¬ 
rium after every plastic instability. The stress-controlled 
path does not equilibrate except when the internal stress 
can balance the external stress. Nevertheless the stress- 
controlled path visits for a while the lowest available 
states until it finds the cheapest escape route (see Fig. E2). 
Of course, increasing the temperature can allow the sys¬ 
tem to escape over higher barriers and thus to open a 
variance with the AQS path. Nevertheless this does 
not change the conclusion that the stress-controlled yield 
is mechanical in nature rather than a thermal melting, 
either real or fictive. This conclusion is strengthened 
by our explicit demonstration that the relative number 
of negative eigenvalues of the Hessian matrix does not 
change as a function of external stress even as a yield 
point is approached. Finally, a clear demonstration that 
the yield point is close to the underlying AQS saddle- 
node bifurcation is provided by the scaling law Eq. HID 
which is found to be observed very accurately in our sim¬ 
ulations. Applications of these ideas to experimental sys¬ 
tem will be reported elsewhere. 
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